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Abstract 



We calculate spin-dependent cross sections for the scattering from mica 
of hypothetical weakly interacting dark-matter particles such as neutralinos. 
The most abundant odd-A isotopes in mica, 27 A1 and 39 K, require different 
shell-model treatments. The calculated cross sections will allow the interpre- 
tation of ongoing experiments that look for tracks due to the interaction of 
dark-matter particles with nuclei in ancient mica. 
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I. INTRODUCTION 



The nature of the dark matter in our galaxy and elsewhere has become increasingly 
puzzling Although it is too early to make definitive statements, ongoing experiments 
seem to imply that there are not enough macroscopic objects in the galactic halo 0] to 
account for the gravitational attraction felt by luminous objects near the edge of the disk. 
One alternative is dark matter in the form of elementary particles. For several years now, 
weakly-interacting massive particles (WIMPs) that arise in supersymmetric extensions of 
the standard model have been attractive candidates @. A variety of experiments to detect 
WIMPs are either already operating or in the planning/prototype stage. 

One promising project Q uses ancient mica as a detector. The idea is that over long 
periods of time a countable number of WIMPs should have collided with underground mica 
and left recognizable tracks. At present, the experiment is not sensitive enough to test the 
WIMP hypothesis fully, but by increasing the amount of mica examined, experimenters hope 
to reach the required level of sensitivity within a few years. Then, confirming or rejecting 
the WIMP hypothesis will require a knowledge of their cross sections for scattering from 
the various elements in mica. In this paper we calculate these cross sections so that the 
experiments can be properly interpreted. Our results are general enough to be applied to 
any heavy Majorana particle; we explicitly calculate the "spin-dependent" response of mica. 
While "spin-independent" cross sections for supersymmetric neutralinos ||, perhaps the 
most plausibly motivated WIMPs, are usually larger than the spin-dependent cross sections 
in most of the nuclei in mica, there are still regions of supersymmetric parameter space 
for which this is not so. Furthermore, WIMPs with no scalar interactions, such as heavy 
Majorana neutrinos, have not yet been completely ruled out. A good calculation of the 
spin-dependent response is therefore warranted. The spin-independent scattering can be 
easily calculated following, for example, the work of Ref . || . 

The composition of mica is 58% 16 0, 16% Si (mostly 28 Si), 12% 27 A1, 5% K (mostly 39 K), 
and 9% H. Since the scattering from even-even nuclei is entirely spin-independent — and 
therefore easily calculable — we need to concentrate here only on 27 A1 and 39 K. (No detailed 
nuclear-structure calculation is required for H.) The nucleus 27 Al is also one of the active 
ingredients in a very high-resolution (AEfwhm ~ 200 eV) and low-threshold [Eth ~ 500 
eV) sapphire-crystal (AI2O3) detector that is currently under development M. The low 
threshold of this experiment makes it ideal for detecting lighter WIMPs. 

Though both aluminum and potassium nuclei are nominally classified to be in the sd shell, 
39 K is at the upper edge of the shell and its wave functions contain significant admixtures 
of 1- and 2- particle-hole states involving higher shells. Most of the correlations in 27 Al, by 
contrast, are produced by purely sd-shell configurations. The techniques we apply in the two 
cases are therefore quite different. After a general discussion of spin-dependent scattering in 
Section |TJ, we will examine the two odd-A nuclei independently in Sections |TT| and [TV]. In 
section [V] we briefly discuss the implications of these calculations for mica-based detection. 

II. SPIN-DEPENDENT SCATTERING 

At low-momentum transfer, the tree-level coupling of neutralinos to quarks (excluding 
the exchange of Higgs bosons) yields the elastic scattering amplitude from a nucleus \N): 
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M = A{N\a p S p + a n S n \N) ■ s x , 



(1) 



where S p and S n are the total proton and neutron spins in the nucleus, s x is the spin of 
the neutralino, a p and a n are coupling constants that depend on the composition of the 
neutralino (in terms of the four or more neutral fermionic partners of the gauge and Higgs 
bosons) and on the distribution of nucleon spin among the different types of quarks, and the 
quantity A depends on other physics at the supersymmetry scale |§. Though these specifics 
apply only to neutralinos, the spin- dependent scattering amplitudes for other Majorana 
WIMPs will have the same general form, Eq. ([I]). The spin-dependent cross section is a 
product of the square of the amplitude and a phase-space factor. 

When the momentum transferred to the nucleus by the WIMP, q, becomes comparable 
to the inverse radius of the nucleus, structure functions (form factors) resembling those 
associated with elastic magnetic electron or neutrino scattering modify Eq. ([!]). The more 
general expression for the spin-dependent cross section is @ 

where Gf is the Fermi weak four-point effective coupling constant, and v the neutralino 
velocity. The subscript A stands for "axial". The function Sa{q) can always be written in 
terms of three fundamental structure functions 

S A (q) = alS 00 (q) + alS n (q) + a aiS 01 (q) . (3) 

The labels and 1 refer to isospin; the constant do is given by a p + a n and at by a p — a n . The 
three structure functions contain expectation values of operators of the form ji(qr) [Ylcx] l±1 
(L even), which depend on spatial coordinates as well as nucleon spins. Details appear in 
Ref. || . Without rewriting all the expressions here, we describe the calculation of Sa(<7) in 
27 Al and 39 K. 



III. ALUMINIUM 

The nucleus 27 Al lies in the middle of the sd shell. To obtain the ground-state wave 
function we diagonalize Wildenthal's "universal" sd interaction || in the full Ohu space. This 
interaction has been meticulously developed and tested over many years and together with 
effective operators accurately reproduces most low-energy observables in sd-shell nuclei |TU 



We carry out our calculations with the Lanczos m-scheme shell-model code CRUNCHER 
IIJ and its auxiliary codes; the m-scheme basis for 27 Al contains 80115 Slater determinants. 
A very similar calculation for 29 Si is reported in |12||, where more details appear. 



The agreement of the calculated spectrum with that measured in 27 A1 is excellent, re- 
flecting the effort that went into constructing the interaction. More relevant is the accuracy 
of the calculated magnetic moment, which derives from the same spin operators that de- 
termine the WIMP structure functions at q 2 = 0. The measured magnetic moment is 
^exp = +3.6415 fix- Using free-particle ^-factors and our wave functions, we obtain a value 
^caic — +3.584 /ijv, which agrees well with experiment. By contrast, the single-particle value 
is Hs P = +4.793 /ijv- 
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Although Ref. |10[ advocates the use of effective operators in the sd shell, we obtain a 
calculated magnetic moment in 27 Al closer to experiment with the free-nucleon spin and 
orbital-angular- momentum operators (the same was true in 29 Si). Unfortunately, we over- 
predict the Gamow- Teller matrix element by a factor of about 1.25, a problem that the 
effective operators would remedy. Thus, we have conflicting evidence on whether to use 
effective operators or not. In the absence of a better prescription, we choose to use the 
free-nucleon g-factors. The consequence for WIMP scattering is that we did not quench our 
calculated cross sections at all. In 39 K, where our calculation provides more information, 
we make a similar decision. 

Ref. []13[ develops two schemes, based on the "Odd Group Model" (OGM) and an ex- 
tension of it (EOGM) for extracting values of neutron and proton spins — the quantities 
that determine spin-dependent WIMP cross sections at q 2 = — from magnetic moments 
and beta-decay lifetimes. These values for 27 Al appear in Table p] alongside the spins and 
orbital angular momenta from our calculation. The most surprising result is the large value 
for (L n ) obtained in the shell-model calculation. On comparing the numbers in Table |, 
we see that the single-particle model (ISPSM) and OGM disagree with the shell model; 
the large value of the neutron's orbital angular momentum is the likely explanation of the 
OGM's failure to reproduce the spin angular momenta. By contrast, the EOGM results 
with gA/dv = 1 (a quenching of the Gamow- Teller matrix element) are very close to our 
shell-model values. This agreement with the most sophisticated phenomenological analysis 
is heartening. However, since spin- dependent neutralino scattering involves the axial nuclear 
current, the failure of any analysis to reproduce beta-decay lifetimes accurately without ad 
hoc quenching — a problem not mentioned in any previous work — is troubling. As a result 
our aluminum calculation has an uncertainty from this source of roughly 30%. The situation 
is better in the potassium calculation since the core-polarization corrections, to be discussed 
in the next section, explain much of the beta-decay quenching. 

As described above, when heavy particles transfer momentum comparable to the in- 
verse nuclear radius, structure functions modify the simple form Eq. ([|) for the scattering 
amplitude. In the shell model, the expectation value of any one-body operator is easy to 
calculate, so the evaluation of structure functions poses no additional problem. We did, 
however, improve upon an approximation in previous works, namely the use of harmonic 
oscillator single-particle wave functions; the true nuclear mean-field potential is closer to a 
Woods-Saxon than a harmonic oscillator form. For comparison we calculate the structure 
functions with both types of wave functions. We use a length parameter of b = 1.73 fm for 



oscillator functions and the standard parameters of Ref. |14[] for the Woods-Saxon potential. 
Fig. |Ij shows the three structure functions in both function of q 2 up to the maxi- 

mum possible momentum transfer for a WIMP moving at 700 km / sec (more than twice the 
mean WIMP velocity). The Woods-Saxon wave functions make only small changes in the 
results. 

For the momentum transfers in the figure, the Woods-Saxon structure functions are given 
to high accuracy by the third-order polynomials 

S 00 = 0.0929516 - 0.472059 y + 1.05996 y 2 - 1.01148 y 3 (4) 
S n = 0.0657232 - 0.449840 y + 1.35041 y 2 - 1.68508 y 3 
Sqi = 0.1563300 - 0.935958 y + 2.45779 y 2 - 2.72621 y 3 , 
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where y = (bq/2) 2 . As noted above, these three functions allow the calculation of spin- 
dependent cross sections for any heavy Majorana particle. 



IV. POTASSIUM 

For the nucleus 39 K, a shell-model diagonalization is difficult to perform because of 
the severe truncations required to the active model space. Since this nucleus is so near 
the boundary between the sd and pf shells, excitations of particles into higher shells can 
have significant effects that are often not well simulated by effective operators. Despite the 
difficulties, our first attempt to treat this nucleus was within the framework of the shell- 
model, with the Hamiltonian and method advocated in Ref . . Computational complexity, 
however, kept us from allowing more than two particles out of the sd shell into the fp 
shell. After diagonalization, the 2p-3h excitations we did include made up approximately 
40% of the ground-state wave function. While the excitations changed the orbital angular 
momentum significantly, they had almost no effect on the spin, and our magnetic moment 
disagreed badly with experiment, indicating that further correlations are indeed important. 
Any further complexity makes use of the shell model difficult, however. We turn instead to 
an alternative scheme, based on perturbation theory, that was successfully implemented a 
few years ago in calculations of several spin- dependent observables in closed-shell-plus (or 
minus)-one nuclei |16| . Since this method has never before been applied to WIMP structure 



functions, we describe it in more detail than we did our sd shell-model results. 

We begin by dividing the Hamiltonian into a one-body term and a residual interaction: 
H = H + V, where H = T + U, the sum of kinetic energy and one-body potential energy 
operators, and V = V — U, the difference between the two- and one-body potential energy 
operators. The eigenfunctions of Ho, which we take to be a harmonic oscillator, form the 
basis of the calculation. The unperturbed ground state of 39 K, a single-proton hole in the d^/2 
orbital, gives the single-particle (Schmidt) magnetic moment: -0.676 /ijv- Corrections come 
from the perturbative expansion in the residual interaction. Very generally, to first order in 
a closed-shell-plus (or minus)-one nucleus the matrix element of a one-body operator, F, is 
given by 

W « HVM+E (M +E «) (5) 

Here a and b are single-hole valence states and a an infinite set of single-hole or two-hole one- 
particle (2h-lp) states constructed from the eigenfunctions of H ; in our work a = b = 0d 3 / 2 - 
The energy denominators also come from Ho; they are integral multiples of 2hu, where %u 
is the characteristic oscillator energy. We choose hu = 45A- 1 / 3 - 25A~ 2 / 3 = 11.1 MeV, 
corresponding to b = 1.934 fm, for A = 39. We explicitly evaluate the intermediate-state 
summations for denominators equal to 2tku, Afwo, 6tku, and 8tku, and then extrapolate the 
results geometrically. 

The two-body matrix elements of V in eq. (Eh contain the one-body potential through the 



—U components, known as Hartree-Fock insertions ||17|| . If the unperturbed Hamiltonian 
Ho were chosen to be the Hartree-Fock Hamiltonian that minimized the energy of a single 
Slater determinant characterizing the closed-shell core, then the Hartree-Fock insertions 
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would have no effect. Ref. [IlBfl shows that even with an oscillator Hamiltonian the Hartree- 



Fock insertions do not affect magnetic moments and their effects ought to be small for 
the operators discussed here. We therefore drop the Hartree-Fock insertions altogether. 
Furthermore, all first-order corrections to magnetic moments, and to the operators used here 
in the q = limit, also vanish because the spin and orbital-angular- momentum operators 
cannot create (or annihilate) a particle-hole state at LS closed shells (selection rules require 
the particle and hole states to have the same orbital structure). At non-zero q, the first-order 
corrections to SA{q), while not identically zero, are still small. It is essential, therefore, to 
consider perturbation theory through to second order, and we add the following additional 
terms to Eq. (||): 

E {ma) MVMWVM E{blFla) («\V\a)(a\V\a) 



(E a — E a )(E a — Efr 



a^a,b 



(E a — E 



\- jffljfflj , , F| > ^ (b\V\b)(b\V\a) , ,,,, , 



,f3^a,b ( E b ~ E fi)( E b ~ E Q 



a^a,b 



(E b - E a y 



(6) 



+ E ^wM-l E { ™iZ lb h ma) 



,^a,b ( E b - E P) 



(E a — E 1 



/3^a,b 



(E a - Ep) 2 



Here, as before, a and b are single-hole valence states, a an infinite set of 2h-lp states, and 
/3, 7 an infinite set of 2h-lp and 3h-2p states. Again the selection rules on (b \ F | a) and 
(a | F | a) forbid the first four terms of eq. @ from contributing at q — 0, and only allow 
small contributions at higher q. The fifth, sixth and seventh terms, however, each contain 
intermediate-state summations over 2h-lp and 3h-2p states that are not constrained by 
selection rules and converge only slowly with increasing energy denominators. [The last two 
terms renormalize the single- hole matrix element (b | F | a) .] These three terms are sometimes 



called the number-conserving set |18j because if the one-body operator were replaced by the 
number operator the terms would contribute nothing. This is apparent when the last three 
terms in Eq. (|6]) are rewritten in the following equivalent form (with a = b): 



Q 
v- 



Q 

F^V 



(7) 



where the operator Q/e is J2b I I /(E a — Ep). If the one-body operator F were to 

commute with both the energy denominator e and the residual interaction V, as does the 
number operator, the second-order correction to the diagonal matrix element would van- 
ish. But spin-dependent one-body operators do not commute with spin-dependent residual 
interactions, resulting, for example, in second-order corrections to magnetic moments from 



tensor forces [19]. 

Here we try two different residual interactions. The first (which we call I), used in 
Ref. [16] for magnetic- moment calculations, is based on the one-boson-exchange potential 
of the Bonn type [ |20|j , but limited to the four or five important meson exchanges. For use 
in finite nuclei, this interaction should be converted to a G-matrix; here we approximate 
the procedure crudely by introducing a short-range correlation function. The resulting 
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interaction has a weak tensor-force component typical of Bonn potentials. Our second 
interaction [21] (called /J here) is a full G- matrix constructed from the Paris potential [22 



and parameterized in terms of sums over Yukawa functions of various ranges and strengths. 
Interaction II has a strong tensor force, producing effects described below. 

Our results for the orbital angular momenta and spins in 39 K appear in Table [TT| As was 
the case earlier, the spins from both forces agree well with those obtained in the phenomeno- 
logical EOGM. To test our results against data, we perform a comprehensive calculation of 
magnetic moments and Gamow- Teller matrix elements, including (under the heading MEC) 
meson-exchange currents, isobar currents, and other relativistic effects (see [[TJJ) in addition 
to the core polarization described above. Table fT\\ presents the results for both isoscalar and 
isovector quantities [The magnetic moment of 39 K is the sum of the isoscalar and isovector 
moments.] The MEC correction is small for isoscalar magnetic moments and beta decay, 
but is significant for isovector electromagnetic operators. There is little difference between 
the two interactions for isoscalar operators and the results are in good agreement with ex- 
perimental values. For the isovector operators, interaction II produces larger corrections and 
poorer agreement with the experimental isovector magnetic moment. But, although neither 
interaction reproduces the Gamow- Teller beta-decay matrix element well, interaction II does 
better, and we recommend it for that reason. Compared with 27 A1, much smaller ad hoc 
quenching is required to reproduce the measured lifetime; the correlations induced by the 
core-polarization produce about two-thirds of the required reduction themselves. 

Turning now to the WIMP structure functions in fig. 2, we see that all three are strongly 
quenched from their single-particle values. The strongest quenching is in Soo(q), which 
is reduced to 25% and 20% of the single-particle value for the two residual interactions 
considered. Our preferred choice, the solid line (corresponding to interaction II), is the 
lowest in all three diagrams. The solid lines are accurately reproduced by the following 
fourth-order polynomials in y = (bq/2) 2 : 

Sqo = 0.0094999 - 0.0619718?/ + 0.162844 y 2 - 0.194282 y 3 + 0.0891054 y 4 (8) 
Sh = 0.0298127 - 0.2176360 y + 0.623646 y 2 - 0.814418 y 3 + 0.4050270 y 4 
S i = 0.0332044 - 0.2319430 y + 0.638528 y 2 - 0.798523 y 3 + 0.3809750 y 4 . 

We repeat that these structure functions completely determine spin-dependent cross sections 
for any heavy Majorana particle. 



V. DISCUSSION 

What are the prospects for mica as a WIMP detector? At present the limits that can be 
set are several orders of magnitude higher than the expected cross sections for neutralinos. 
But by analyzing more mica, the experiments may reach the necessary level of sensitivity 
within a few years. Because the efficiency with which each element in mica can be detected 
when it recoils is different, we cannot easily present an overall cross section for a representa- 
tive neutralino, as is the common practice in the literature. But we note that the calculated 
efficiencies [E31 begin to increase at a recoil energy of about 20 keV, corresponding to a q 2 



of about 10" 3 GeV 2 in Al, and about 1.6 x 10" 3 GeV 2 in K. If the 20 keV were a detection 
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threshold, then as the curves in figs. 1 and 2 show, the number of events below threshold 
would be a small fraction of the total. 

Whatever the precise characteristics of the detector, the expected number of events 
can now be calculated by folding the efficiencies together with the cross sections presented 
here and with the flux of WIMPs, which, unfortunately, is still poorly constrained. But 
continuing work with the mica, with complementary detectors, and with telescopes should 
make it possible to rule in or out most hypothetical WIMPs within the next decade. 
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FIGURES 



FIG. 1. The three 27 A1 structure functions Soo (top), Sn (middle), and Soi (bottom) as a 
function of momentum-transfer squared. The dotted lines are the single-particle functions, the 
dashed lines the full functions with harmonic oscillator single-particle wave functions, and the 
solid lines the full functions with Woods-Saxon single-particle wave functions. 

FIG. 2. The three 39 K structure functions Soo (top), Sn (middle), and Soi (bottom) as a 
function of momentum-transfer squared. The dotted lines are the single-particle functions, the 
dashed lines the full functions obtained from the modified Bonn interaction, and the solid lines the 
full functions obtained from the Paris-based G-matrix (see text). 
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TABLES 



TABLE I. The value of (the z-projection of) the nuclear spin and orbital-angular- momentum 
matrix elements for 27 A1. 





(S P ) 




(S n ) 


(Lp) 


<L„) 


ISPSM 


0.5 







2.0 





OGM 


0.25 







2.25 





EOGM (g A /g v = 1.0) 


0.333 




0.043 






EOGM (g A /g v = 1.25) 


0.304 




0.073 






Shell Model 


0.3430 




0.0296 


1.7812 


0.3461 


TABLE II. The value of (the z-projection of) the nuclear spin and orbital-angular-momentum 


matrix elements for 39 K. 














(S P ) 




(S n ) 


(L P ) 


(L n ) 


ISPSM 


-0.3 







1.8 





OGM 














EOGM (g A /gv = 1.0) 


-0.196 




0.055 






EOGM {g A /gv = L25) 


-0.171 




0.030 






Force I 


-0.197 




0.051 


1.339 


0.308 


Force II 


-0.184 




0.054 


1.068 


0.564 


TABLE III. Corrections to the single-particle (sp) values of the isoscalar (/j,g) 


and isovector 


(/xy) magnetic moments and diagonal Gamow- Teller beta-decay matrix element (M(GT)), for the 


mass A = 39 system from a 


Bonn-based, weak-tensor force, I, and 


a Paris-based, strong-tensor 


force, II. The corrections 


are due primarily to second-order core- 


-polarization (CP) and me- 


son-exchange currents (MEC). See text for further 


explanation. 






Force 


sp 


CP 


MEC 


Sum 


Expt 


US I 


0.636 


0.066 


0.004 


0.706 


0.707 


II 


0.636 


0.064 


0.008 


0.708 


0.707 


Hv I 


-0.512 


-0.138 


0.364 


-0.286 


-0.315 


II 


-0.512 


-0.378 


0.369 


-0.521 


-0.315 


M(GT) I 


-0.976 


0.176 


-0.009 


-0.809 


-0.647 


II 


-0.976 


0.205 


-0.016 


-0.787 


-0.647 



11 



